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Abstract 

We study the charged t — ► bWh° and t — ► bWA° decays in the framework of the 
general two Higgs doublet model, so called model III and beyond. Here, we take the 
Yukawa couplings complex and introduce a new complex parameter due to the physics 
beyond the model III, to switch on the CP violating effects. We predict the branching 
ratios as BR(t bWh°) ~ 1CT 6 and BR(t bWA°) ~ 10^ 8 . Furthermore, we observe 
a measurable CP asymmetry, at the order of 10 -2 , for both decays. 
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1 Introduction 



Because of its large mass, the top quark has rich decay products and this opens a new window 
to test the standard model (SM) and to get some clues about the new physics, beyond. In 
the literature there are various studies in the SM and beyond []I||-||13|.The rare flavor changing 
transitions t — > cg(y,Z) have been studied in [ffl H], t — > cH° in §], 0, |§, |9| and t — ► 
c/1/2 in ||10|| . The SM predictions of the branching ratio (BR) of the process t — > cg(y,Z) 
is 4 x 1CT 11 (5 x 1CT 13 , 1.3 x 1CT 13 ) [@, and t — > cH° is at the order of the magnitude of 
1CT 14 — 1CT 13 , in the SM 0, which are not measurable quantities even at the highest luminosity 
accelerators. Possible new physics effects are the candidates for the enhancement of the BR's 
of the above processes, t — > cH° and t — > c/1/2 decays have been analysed in |9j and [Hj], 
in the framework of the general two Higgs doublet model (model III). In these studies it has 
been observed that there could be a strong enhancement in the BR, almost seven orders larger 
compared to the one in the SM for the decay t — > cH°; a measurable BR, at the order of the 
magnitude of 10~ 8 — 1CT 7 for the decay t — > 0/1/2- in [01 t — > cVV decay has been analysed in 
the topcolor assisted technicolor theory. 

The charged t — > b transitions exist in the SM model and have been studied in the literature 
extensively. The top decay t — ► bW has been analysed (see |T2| and references therein) in the 



two Higgs doublet model and t — > bWZ decay has been studied in p3| . 

The present work is devoted to the analysis of the charged t — > bW h° and t —>■ bW A 
decays in the framework of the general two Higgs doublet model (model III). This decay occurs 
in the tree level with the extended Higgs sector since the scalar bosons h° and A exist in the 
new sector. We study the BR of the above decays and obtain a measurable quantities, at the 
order of the magnitude of 10~ 6 and 10~ 8 , respectively. Furthermore, we search the possible 
CP violating effects. To obtain a nonzero CP asymmetry Aqp we take Yukawa coupling for 
bh° (A ) b transition complex and introduce a new complex parameter, where its complexity 
comes from some type of radiative corrections, due to the model beyond the model III (see 
section II). We obtain a measurable Aqp, at the order of the magnitude of 10~ 2 and observe 
that these physical quantities can give valuable information about physics beyond the SM, the 
free parameters existing in these models. 

The paper is organized as follows: In Section 2, we present the BR and Aqp of the decay 
t — > bW h°(A°) in the framework of model III. Section 3 is devoted to discussion and our 
conclusions. 
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2 t — > and t — > bWA° decays with possible CP vio- 

lating effects. 

If one respects the current mass values of h°(A°), namely m h o ~ 85GeV (m^o ~ 90GeV), the 
charged t — ► 6W / /t°(yl ) is kinematically possible and does not exist in the SM model. With the 
minimal extension of the Higgs sector the CP odd new Higgs scalar A arises and the t — ► bWA° 
decay in the tree level is permitted. In this model, t — > bWh° decay is possible in the tree level, 
where h° is the new CP even Higgs scalar and, in general, it mixes with the SM one, H°. In 
this section, we study the BR in the general two Higgs doublet model, so called model III and 
the possible CP violating asymmetry, beyond. 

The t — > bWh°(A°) decay is created by the charged t — > bW process and the neutral 
t — > t*h°(A°) or b* — > bh°(A°) processes, which are controlled by the Yukawa interaction 

Cy = vfjQiJiUjR + nl'QnOiPjH + C',j Q,,.o,r in + ^QafcDjR + h.c. , (1) 

where L and R denote chiral projections L(R) = 1/2(1 =f 75), (pi for i — 1,2, are the two scalar 
doublets, QiL are left handed quark doublets, Ujn(Dj R ) are right handed up (down) quark 
singlets, with family indices The Yukawa matrices r^- D and £^- D have in general complex 
entries. By considering the gauge and CP invariant Higgs potential which spontaneously breaks 
SU{2) x U{1) down to 17(1) as 

VfaM = c 1 (0+0 1 -^ 2 /2) 2 + c 2 (0 2 + 2 ) 2 

+ - V 2 /2) + 0+0 2 ] 2 + CiMMtffa) ~ (0r0 2 )(0 2 + 0l)] 

+ c 5 [i?e(0+0 2 )] 2 + c 6 [7m(0+0 2 )] 2 + c 7 . (2) 



and choosing the parametrization for <f>i and 2 as 

\, ( V2 X + " 



01 = 71 



;02 = ^(^^J- (3) 



with the vacuum expectation values, 



<0i>=7^(°) ;<0 2 >=o, (4) 

the Hi and H 2 becomes the mass eigenstates h° and A respectively since no mixing occurs 
between two CP-even neutral bosons H° and h°, in tree level. This scenerio permits one to 
collect SM particles in the first doublet and new particles in the second one. Furthermore the 
Flavor Changing (FC) interaction can be obtained as 

£y,fc = ^QiL&UjR + ^QiLfoDjR + h.c. , (5) 
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with the couplings £ U,D for the FC charged interactions 

£ch = £n v ckm , 

ich — Vckm £n > (6) 
where is defined by the expression 

er } = {v^re^v^ . (7) 

Notice that the index "N" in t^ D denotes the word "neutral". 

Using the relevant diagrams for the t — > bWh°(A°) decay which are given in Fig [l] and taking 
into account only the real Yukawa couplings £j},bb> £n,ui the matrix element square \M\ 2 {h°) 
{\M\ 2 {A )) reads 

\M\lo {AO) (pl,p b , k, q) = A(h° (A )) + (Z° bb f f 2 (h° (A )) + (^ tt ) 2 f 3 (h° (A )) (8) 

where 

h(h Q ) = lQ\V tb \ 2 m b m t (m 2 w (s 2 2 (h )-s 2 1 (h ) + 2s 1 (h°)s^^ 
+ 2 Sl (h°) (( Sl (h°) - s 2 (h )) k.{p 1 -p b ) + 2s 1 {h°)p 1 .p b ) 
+ 4- f - ( s l(h°) + 2 sl{h°) + 2 Sl (h°) s 2 {h )) (k.q) 2 + 2 8l {h°) (2 Sl {h°) 

+ s 2 (h )) k.qq.(p 1 - p b ) + 8 s\(h Q ) p x .q p b .q)) , 

Mh°) = 8\V tb \ 2 ^-m 2 w (s 2 (h°) + Sl (h c 

+ J_ k .q((s 2 (h )+2 Sl (h )) 

x (s 2 (h°) k.qpx.p b + 2si(/i°) k.p b q.pt) -2si(/i°) s 2 (/i°) k.pip b .qj 
+ 2 (Ai.pi A;.p 6 - z h o q.pi q.ph)^ , 

/ 3 (^°) = 8 \V tb \ 2 ( - m 2 w (A Sl (h°) ( Sl (h°) - s 2 (h )) x t k. Pb + (( Sl (h°) + s 2 (h )) 2 x h o 
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X h opx.p b 



- 4sl(ti } )x t )pi.pb 

+ 4r k.q (s 2 (h°) (s 2 (h°) + 2 Sl (h )) k.q Pl .p b - 2 8l (h°) s 2 (h°) q. Pl k.p b + 2 Sl (h°) 
x (2s 1 (h°)+s 2 (h°))k.p 1 q.p b 

+ 2 Sl (h°) ( Sl (h°) k. Pl k.p b - (2 (2 Sl (h°) + s 2 (h )) x t k.q + Sl (h°) (x h o - Ax t ) q.p, ) 



= 16 \V tb \ 2 m b m t (m 2 w (s 2 2 (A°) + sl(A )) x A o + -L-(2^(A°) - sj(A )) (k.q) 2 ^j , 



f 2 (A°) = 8\V tb \ 2 ^-m 2 w (s 2 (A°) + s 1 (A )) 2 x A op 1 . Pb 

+ \ k.q ((s 2 (A°) + 2 Sl (A )) (s 2 (A°) k.q Pl .p b + 2 Sl (A°) k. PbPl .q) 

- 2s 1 (A°)s 2 (A°)k.p 1 p b .q) +2s 2 (A°)(k.p 1 k.p b -x A oq.p 1 q.p b ) ) J , 

f 3 (A°) = 8\V tb \ 2 ^-m 2 w ( Sl (h )-s 2 (h )yx A o Pl .p b 

+ 4" k -Q ( S ^ A °) ~ 2 s ^ h °)) k-QPi-Pb + 2 s 1 (A°) (2 Sl (A°) - s 2 (A )) k. Pl p h .q) 

+ 2s l (A°)s 2 (A°)k.p b p 1 .q) +2s 2 l (A°)(k.p 1 k.p b -x A op l .qp b .q)y (9) 

Here the functions si( 2 ,3)(/i (A )) are 



si(h°) 



9w 



4m 2 w (l+x t -2^)> 
9w 



S2{k0) 2/,/f r (l • x h o -y t 2 '■ 



sm(A°) = (-) s m (h° -+ A°) , (10) 

with weak coupling constant g w , x h o( A o) = ^ A ° } , x t = ^f- and y t = -^r- and pi, p?,, g and k 
are four momentum of t quark, b quark, W boson and Higgs scalar h°(A°), respectively. 
Finally, using the well known expression defined in the t quark rest frame 

(2tt) 4 (4) d 3 p b d 3 q d 3 k 

dT h o {AO) = — b K >{pi - Pb - k - q) 



12 mt v ' (2n) 3 2E b (2ir) 3 2E w (2vr) 3 2 E h o (AO) 

x \ M \lo(AO)(Pi>Pb,k,q) (11) 

and the total decay width T T ~ T(t — > 514 7 ) as T T = 1.55 GeV, we get the BR for the decay 
t -> 

Now, we would like to study a possible CP violating effects, which can give comprehensive 
information about the free parameters of the model used. For the process under consideration, 
the CP violation can be obtained by choosing the complex Yukawa couplings in general, namely, 
taking the parametrizations 

fU __ \pV I iOtt 

^nm = \€nm\ e%6tb ■ (12) 
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However, this choice is not enough to get non-zero Acp 

r - f 

Acp = (13) 

where f is the decay width for the CP conjugate process. This forces one to go beyond the 
model III and try to obtain a new complex quantity so that its complexity does not come from 
the Yukawa couplings but from some radiative corrections. Under the light of this discussion, 
we introduce an additional complex correction x to b — > b transition, which may come from the 
new model beyond the model III as 

(£,N,bb + ^N,bb) + (^N,bb ~ ^N,bb)l5 + X 

and 

(£,N,bb ~ £.N,bb) ~ (Cv,M> + Cv,t*)75 + X75 > 

Here we take the magnitude of x & t most \x\ ~ 1CT 2 , which is more than one order smaller 
compared to the vertex due to model III. In this case, we take the correction to the t —*■ t 
transition small since the strength of t — > t transition is weaker compared to strength of the 
b — > b transition, with respect to our choice (see Discussion section). 

At this stage, we introduce a model beyond the model III as follows: The multi Higgs 
doublet model which contains more than two Higgs doublets in the Higgs sector can be one 
of the candidate. The choice of three Higgs doublets brings new Yukawa couplings which are 
responsible with the interactions between new Higgs particles and the fermions. The Yukawa 
lagrangian in three Higgs doublet model (3HDM) reads 

+ p^QiLfoDjR + h.c. , (14) 
where p^ D ^ is the new coupling and 03 can be chosen as 



J_ / V2F+ 



(15) 



with vanishing vacuum expectation value. The fields F + and H s (H4) represent the new charged 
and CP even (odd) Higgs particles, respectively. Notice that the other Yukawa couplings 
and Higgs particles in eq. (|T4| ) are the ones existing in the model III. Now, we choose the 
additional Yukawa couplings p^ D ' ) real and take into account the radiative corrections to the 
b — > b transition which comes from the contributions of third Higgs doublet for the decay 
under consideration. Here the complexity of the parameter should come from the radiative 



5 



corrections but not from the new Yukawa couplings. We can take this complex contribution as 
a source for the additional part x- Since the number of free parameters, namely masses of new 
Higgs particles m F ±, m H:i , m Hi and the new Yukawa couplings p^ D \ increases, there arises a 
difficulty to restrict them. However, the overall uncertainity coming from these free parameters 
lies in the contribution x an d it can be overcome by the possible future measurement of the 
CP violation for our process. 

Finally, by using the definition 

d 2 r(t^bWh°(A )) _ d 2 r(t^bWh°(A°) 
A h°(A°) (F F \ dE b dE w dE b dE w 

^cp \.^w,-^b) d 2 r(t^bWh°(A°)) | d?r(t-+bWhO(AQ) ^ 

dE b dE\y dE b dEw 



we obtain the differential Acp(Ew, E b ) for the process t — > bWh°(A°) as 

Acp A °\E w , E b ) = \^ bb \ \x\ sin6 bb st^ ^^ , 



(17) 



where 
$ fc0 = 



4m t Sl (h°) \V tb \ 2 (4(2 E 2 W m t si(/i°) (x t - 2) - E 2 (2 E w + m t ) s 2 {h°) (1 + x t ) 

+ E b E w (E w (s 2 (h°)(l + 3x h o -3x t )+4 Sl (h )(l + 2x h o)) 

+ m t (s 2 (h ) + 2s 1 (h°)(2x h0 +x t )))) 

+ m 2 w (m t (s 2 (h°) (-1 + (x h o - x t ) 2 ) + 4 Sl (h°) (l + x t + x h o) + 2 E b (-2 Sl {h°) (1 + 2 x h o + x t ) 

+ s 2 (h°)(-l + x h o -x t )(2x t -l)) 

- 2 E w (s 2 (h°) (x h o - x t ) (l-x t + x h o) + Sl (h°) (4 + 2 x h o (2 + 2 x h o -x t ) + 2 x t (3 - x t )))) 

+ Jr (E b E 2 w (E b (2E w + m t )s 2 (h Q ) - A E w m t Sl {h°)))j , 

$ A ° = 4m tSl (A°) \V tb \ 2 (^{2E 2 w m tSl (A ) (x t - 2) - E 2 (2 E w + m t ) s 2 (A°) (1 + x t ) 

+ -Eft -EW 

(E w {s 2 {A°){\ + 3x A o - 3x t ) +4si(A°)(l + 2x A0 )) 

+ m t ( S2 (A°) + 2 Sl (A°)(2x A o+x t )))) 

+ m^(-m t (s 2 (A°) (-1 + (x A o -x t ) 2 )+A Sl A° (l + x t + x A o) 

+ 2 £ 6 (2 si(A°) (1 + 2 x A o + x t ) + s 2 (A°) (-1 + x t + 2 x\ + x A o(l - 2 x t )) 

+ 2E W (s 2 (A°) (x A o - x t ) (l-x t + x A o) + Sl (A°) (4 + 4 x A o (1 + X40 - 2 x t ) + 2 x t (3 - x t )))) 

+ -^(E b ^( J E 6 (2^ + m t ) S2 (A )-4Ep F m tSl (A )))) , (18) 



with x — e 10 * \x\i ^N,bb = e%dhh \^N,bb\- Notice that we do not present the functions D(h°) and 
D(A°) since their explicit expressions are long. Here the functions si(h°(A )) and S2{h°(A )) 
are given in eq.fllQl). 



3 Discussion 

This section is devoted to the analysis of the BR and Acp of the decay t — > bWh° and t —>■ bWA° 
in the framework of model III and beyond. In our numerical analysis we use the form of the 



coupling f$r%\ which is defined as = yJ^fC^ij ■ 

Since the model III contains large number of free parameters such as Yukawa couplings, 
%Nij i ^ e masses °f new Higgs bosons, H ± , h° and A , we try to restrict them by using 
experimental measurements. In our calculations, we neglect all the Yukawa couplings except 
£jy tt and due to their their light flavor contents. In addition to this we neglect the off 

diagonal coupling £jy te , since it is smaller compared to £j^ ti (see fll4l). One of the most important 



experimental measurement for the prediction of the constraint region for the couplings £jy tt and 
^Nbb * s the the CLEO measurement |15| 

BR(B -> X s7 ) = (3.15 ± 0.35 ± 0.32) 10~ 4 . (19) 

and our procedure is to restrict the Wilson coefficient Cj which is the effective coefficient of 
the operator 7 = -^^s a a flv {m h R + m s L)b a T^ v (see JTJJ and references therein), in the region 



0.257 < \Cj ff \ < 0.439, where the upper and lower limits were calculated using eq. (|19f) and 
all possible uncertainities in the calculation of Cj^ [f[4|| . In the case of the calculation of Aqp, 
£.Nbb (Cnu) i s taken complex (real) and a new small complex parameter x, due to the physics 
beyond model III, is introduced. In the following, we choose \r t b\ = \ jrr^\ < 1- Notice that, in 
figures, the BR and Aqp are restricted in the region between solid (dashed) lines for > 
(C?^ < 0). Here, there are two possible solutions for Cj due to the cases where |r#,| < 1 and 
Tth > 1. In the case of complex Yukawa couplings, only the solutions obeying \r t b\ < 1 exist. 

In Fig. @, we plot the BR(t -> bWh°) with respect to ^ for m H ± = 400 GeV, m h o = 
85 GeV . As shown in this figure, the BR is at the order of the magnitude of 10~ 6 and it 
increases with the increasing values of the Its magnitude (the restriction region ) is larger 
(broader) for Cj* > compared to the one for Cy' < 0. 

Fig. |3] is devoted to the same dependence of the BR(t — > bWA°) for m H ± = 400 GeV, 
m A° — 90 GeV. For this process the BR is at the order of the magnitude of 10~ 8 , almost 2 
order smaller compared to the BR(t —>■ bWh°). It increases with the increasing values of the 
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and its magnitude (the restriction region) is larger (broader) for Cy > compared to 

f f f° 
the one for Cj < 0. Furthermore, the restriction region is sensitive to the parameter and 

for C e 7 ff < 0, upper and lower bounds almost coincide. 

Fig. H (||) represents BR(t — > 6H^/t°(A )) with respect to m h o{m A a) for m#± = 400 GeV^ 
and £jv 6b = 30m&. Here the 5i? increases with the decreasing values of 771^0(771^0). This can 
give a powerfull information about the lower limit of the mass value 777^0(771^0) with the help 
of the possible future experimental measurement of the process under consideration. Notice 
that with the increasing values of 777^0(777^0) the restriction regions for C 7 e// > and C 7 e// < 
become narrower and coincide. 

Now, we would like to analyse the CP asymmetry Aqp of the decay t — * bWh°(A°). To ob- 
tain a nonzero Aqp we take the coupling ^ bb complex and introduce a new complex parameter 
X due to the physics beyond the model III (see section II). 

In Fig I (|7|) we present the sinO^ dependence of Acpit — > bWh°(A )) for \^Nbb\ = 30 m&, 
|x| = 10~ 2 , the intermediate value oisin6 x = 0.5 and m h o = 85 GeV {m A o = 90 GeV). A C p is at 
the order of the magnitude of 10 -3 — 10 -2 and slightly larger for Cj^ < compared to the one for 
C e 7 ff > 0. Fig |(|) represents the sin0 x dependence of A CP (t -> bWh°(A )) for = 30 m b , 

\x\ = 10~ 2 , the intermediate value of sin9bb = 0.5 and m h o = 85GeV(mA0 = 90 GeV). The 
behavior of Aqp is similar to the one obtained in Fig. ^| (0). As shown in these figures, Acp is 
a mesurable quantity, which gives strong clues about the possible physics beyond the SM. 

Now we will summarize our results: 

• The BR of the process t — > bWh°(A°) is at the order of 10~ 6 (10~ 8 ) in the model III and 
it can be measured in the future experiments. This ensures a crucial test for the new 
physics beyond the SM. 

• The BR is sensitive to £jv(,& and the mass value m^o (m^o). This is important in the 
prediction of the lower limit of the mass m h o(m A o) with the possible future experimental 
measurement of the process under consideration. 

• Aqp is at the order of the magnitude of 10~ 2 for the intermeditate values of sin O^b and and 
sin9 x . This a mesurable quantity, which gives a strong clue about the possible physics 
beyond the SM and also more beyond. 

Therefore, the experimental investigation of the process t — > bWh°(A°) will be effective for 
understanding the physics beyond the SM. 
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( c ) 



Figure 1: The diagrams contribute to the decay t — > bWh°(A°). 
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15 20 25 30 35 40 



Figure 2: BR(t -> bWh°) as a function of 2£» for m^± = 400 GW, m ft0 = 85 Gel/. Here 

the is restricted in the region bounded by solid lines for Cy > and by dashed lines for 
C 7 e// <0. 
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Figure 3: The same as Fig. |]but for the decay BR{t — > bWA°). 
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Figure 4: BR{t — > bWh°) as a function of m h o for ^ bb = 30 mf,, m^± = 400 GeV. Here the 

BR is restricted in the region bounded by solid lines for > and by dashed lines for 

C-/ J ■ 0. 
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Figure 5: The same as Fig. f|but for the decay BR(t — > bWA°). 
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sin 8bb 



Figure 6: A C p(t — > bWh°) as a function of sin9 bb for \^ bb \ = 30 m b , m H ± = 400 GeV, 
\x\ = 10~ 2 , sin# x = 0.5. Here the A C p is restricted in the region bounded by solid lines for 
Cj > and by dashed lines for < 0. 




sin 9bb 



Figure 7: The same as Fig. | but for the decay A CP (t -> 6WM°). 
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sin6 x 



Figure 8: A C p{t — > bWh°) as a function of sin9 x for |^ 66 | = 30 m&, m#± = 400 GeV^, |x| 
10 -2 , sin9bb = 0.5. Here the Acp is restricted in the region bounded by solid lines for > 
and by dashed lines for < 0. 
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Figure 9: The same as Fig. ^| but for the decay A C p(t — > 
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